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t-h Abstract 
O 

We calculate the md-m u contribution to the mass splittings in baryonic isospin multiplets using 
SU(3) chiral perturbation theory and lattice QCD. Fitting isospin- averaged perturbation theory 
functions to PACS-CS and QCDSF-UKQCD Collaboration lattice simulations of octet baryon 
masses, and using the physical light quark mass ratio m Ll /m d as input, allows M n - M p , M z - M z + 
and Ms- - M^o to be evaluated from the full SU(3) theory. The resulting values for each mass 
splitting are consistent with the experimental values after allowing for electromagnetic corrections. 
In the case of the nucleon, we find M n -M„ = 2.9+0.4 MeV, with the dominant uncertainty arising 
i from the error in m u /m d . 
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00 1- Introduction 

The physical mass splittings between members of baryonic isospin multiplets have been mea- 
sured extremely precisely flUEl: 



x 



M„ -M p = 1.2933322(4) MeV, (la) 
M z - - M z+ = 8.079(76) MeV, (lb) 
M E - M E o = 6.85(21) MeV. (lc) 

However, the decomposition of each into its two components, arising from electromagnetic effects 
and the d - u quark mass difference, is less well known. Clearly, once one contribution has been 
well determined the other can be inferred from the total. 

In recent years, several research groups have presented lattice determinations of both the 
QCD contribution to the baryon mass splittings, e.g., J3]-|6|, and the electromagnetic contribu- 
tion, e.g., [T7T[T0l. This work uses SU(3) chiral perturbation theory to determine estimates of the 
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Table 1: Coefficients of terms linear in the non-strange quark mass, Bmi — > and the strange quark mass, 

6m, — > (m^. - ml/2), expressed in terms of the leading quark-mass insertion parameters a, /? and cr. 



strong contribution, by fitting currently available isospin-averaged lattice calculations IfTTl IT21 of 
the octet baryon masses. We do not consider electromagnetic effects. 

Our evaluation of the strong contribution to M n - M p is of particular interest in the light of 
recent results which suggest that the accepted value for the electromagnetic contribution, namely 
A EM = 0.76 + 0.30 MeV 03, may be too small. Walker-Loud et al. (WLCM) claim to find 
an omission in the traditional analysis and present a larger value of 1.3 + 0.6 MeV 031- From 
these estimates, one infers strong isospin breaking contributions of A mrf _„, i( = 2.05 + 0.30 MeV 
(traditional) and 2.6 + 0.6 MeV (WLCM) respectively. Clearly, independent theoretical estimates 
of the size of the strong contribution to M n - M p , such as reported here, are of considerable value. 

2. Method 

2.1. Isospin-averaged fit 

The fit to isospin-averaged PACS-CS lattice results p"T|| which we use for this work has been 
reported on in previous papers; we refer to 0314T71I for details. In brief, we use a standard heavy- 
baryon chiral perturbation theory formulation with a finite-range regularization scheme (FRR), 
discussed in [fT8l - l22ll . The mass of a baryon B in the chiral expansion is written as 

M B = M (0) + 5Mf + 5Mf 2) + (2) 

where the leading term M (0) denotes the degenerate mass of the octet baryons in the SU(3) chiral 
limit and is independent of the quark mass matrix M q and baryon B. The notation denotes 
the contribution to the octet baryon mass at order M^'\ 

The correction linear in the quark masses can be expressed as 

SMf = -C^Bmi - C ( B l) s Bm s , (3) 

with the coefficients given in Table [IJ The leading loop corrections which contribute to 8M^ IT> are 
made explicit in [16] and include both octet and decuplet baryon intermediate states. 

We retain the octet-decuplet mass splitting 8 in numerical evaluations, setting 6 = 0.292 GeV 
to the physical N - A splitting. The baryon-baryon-meson coupling constants are taken from 
phenomenology; D + F = g A = 1.27, F = \D and C = -2D, and / is set to / = 0.0871 GeV, 
a chiral perturbation theory estimate for the pion decay constant in the SU(3) chiral limit Il2"3"l . 
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The fit parameters are the octet baryon mass in the chiral limit M (0) , the SU(3) chiral symmetry 
breaking parameters a, 8, cr, and the finite-range regulator mass A, and correspond to those in 

Refs. nana. 

The fit to the PACS-CS baryon octet data is shown in IfTTll . and a comparison between the 
experimental values and the octet baryon masses evaluated at the physical point is given in Table[2} 
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0.959(24)(9) 


0.939 
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1.129(15)(6) 


1.116 
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1.188(H)(6) 


1.193 




1.325(6)(2) 


1.318 



Table 2: Extracted masses for the octet baryons. The first uncertainty quoted is statistical, while the second allows 
for variation of the form of the UV regulator and 10% deviation of /, F, C, and 6 from their central values. The 
experimental baryon masses are shown for comparison. 



2.2. Evaluation ofSU(2) mass splittings 

A feature of SU(3) chiral perturbation theory is that the same coefficients appear in the baryon 
mass expansion both when including and excluding SU(2) breaking effects. These coefficients, 
once evaluated by fitting to the i so spin- averaged lattice results as done above, can thus be used 
to provide information about the SU(2) breaking mass splittings. That is, we can use Nt = 2 + 1 
flavor lattice simulations, which are currently available, and are computationally cheaper than the 
SU(2)-broken Nf = 1 + 1 + 1, to derive our results. 

To calculate the baryon mass splittings, we modify the SU(3) chiral perturbation theory expan- 
sions used in the previous section to allow for a non-zero light quark mass splitting: m d - m u £ 0. 
While this is a straightforward extension, we note that it generates a n°T] mixing term in the SU(3) 
Lagrangian. The fields must thus be diagonalized into the mass basis via a field rotation. 

To be explicit, we recall the usual definition of the meson field: 
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and meson Lagrangian: 



f 2 

Xeff = yTrO^Z 1 "^) + ATr(M q (J? + £)), 



(5) 



where M q = diag(m M , m^, m s ) is, as above, the quark mass matrix. 
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Expanding this Lagrangian in powers of the meson field, the mass term can be written as 

X kin =5Tr(M„0 2 ) (6a) 



=B(m u + m d ){n + n-) + B(m s + m d ){K v K) (6b) 
+ B(m s + m u )(K + K-) + ^(m u + m u )((7r ) 2 ) 

+ f (m d + m u + Am s ){r] 2 ) + -^=(m u - m d )(7]n°). 
6 V3 

Clearly, for m u ^ m d , mixing occurs between the 7r° and 77. 

To identify the meson masses we remove this mixing and bring the kinetic term into the canon- 
ical form via a field rotation: 

7T° — » 7T°COS6 - 77sin6, (7a) 
r] —¥ 7r°sine + 77COS6, (7b) 



where the mixing angle e is given by 



tan2e = — ^ - — (8) 
2m s - (m d + m u ) ' 



After performing this rotation, the SU(3) meson masses take the form: 
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(9e) 



where m ff o and m n now contain some dependence on the mixing angle e. 

This extension generates a separate mass expansion, of the form of Equation [2} for each mem- 
ber of the baryon octet. The terms linear in quark mass can be expressed as 

6Mf = -Cf u Bm u - C { l\Bm d - C%Bm s , (10) 

where the coefficients Cs q are given explicitly in Table [3} 

The loop contributions SM^ /2 ' have the same form as in the i so spin- averaged case, with sep- 
arate couplings and integrals for each of the mesons it* ,n° ,K ± ,K° ',?]. The tt* and K* remain 
pairwise mass-degenerate. Of course, because of the redefinition of the meson fields, the baryon- 
baryon-meson couplings will also receive contributions depending on e, and are now complicated 
functions of quark mass and the coupling constants F, D, and C. These are given in Table s|4] and [5} 
As expected, setting e — > returns the usual isospin-averaged functions. 
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Table 3: Values for the terms linear in the up, down and strange quark masses, expressed in terms of the SU(3) 
breaking parameters a, ft and cr. 
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Table 4: Chiral SU(3) coefficients for the coupling of the octet baryons to decuplet (T) baryons through the pseu- 
doscalar octet meson (p. 



It is now straightforward to write expressions for the baryon mass splittings as a function of 
quark mass only. All other free parameters, namely the SU(3) breaking parameters a, B and <x, as 
well as the regulator mass A, are specified by the isospin-averaged fit described previously. 

To evaluate the baryon mass splittings at the physical point, we input the physical light-quark 
mass ratio R := — . The Gell-Mann-Oakes Renner relation suggests the definition 



B(m d -m u ) 1(1 -R) 



u = v """„ : = ~ 77 . Tz m* , (11) 



2(1+1?) 



*(phy S ) ' 



which allows us to define 



Bm u = m n(vbys) /2-(0, (12a) 
Bm d = m 2 (phys) /2- u>, (12b) 
Bm s = ml -ml /2. (12c) 

4 "(phys) ^(phys)' V ' 

Here, we take w„- (phys) = 137.3 MeV and m K(phys) = 497.5 MeV to be the physical isospin-averaged 
meson masses [0. 

Evaluating the mass splitting expressions at these 'physical' quark masses, with loop meson 
masses calculated using Equation |9} then gives our estimate of the baryon mass differences at the 



physical point. A discussion of the error analysis is given in section 3.1 
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Table 5: Chiral SU(3) coefficients for the coupling of the octet baryons to octet (B) baryons through the pseudoscalar 
octet meson (p. 



3. Results 

The calculation outlined in the previous section, where the central value is obtained using the 
stated phenomenological estimates for /, F, C and 6, gives 

M n -M p = (w/m 2 (phys) )(20.3 + 1 .2) MeV (13a) 
M x - - M z + = (a;/m 2 (phys) X52.6 ± 2.0) MeV (13b) 
M 3 - - M E o = (o;/m 2 (phys) )(32.3 + 1.6) MeV. (13c) 

The quoted uncertainties contain all statistical and systematic errors, discussed in the following 
section, combined in quadrature. 

We take as input the current best-estimate for the physical up-down quark mass ratio Il2"4"ll . 

in 

R= -JL = 0.553 + 0.043, (14) 

m d 



determined by a fit to meson decay rates. We note that this value is compatible with more recent 
estimates of the ratio from 2 + 1 and 3 flavor QCD and QED 0, [251 • Including in quadrature the 
uncertainty due to the stated range for R, we find 



M n -M i 
- M s 



' p = 2.9 ± 0.4 MeV 
:+ = 7.5 + 1.0 MeV 



(15a) 
(15b) 
(15c) 



M s - - M^o 



= 4.6 + 0.6 MeV. 



3.1. Statistical and systematic uncertainties 

The errors quoted are the result of a complete error analysis, taking into account the correlated 
uncertainties arising from all of the fit parameters, as well as propagating the quoted uncertainty 
in R. We estimate the systematic error in our result by considering variations of the regulator and 
allowing for deviation of the phenomenologically set parameters /, F, C and 6 from their central 
values by +10%. 

Monopole, dipole, Gaussian and sharp cutoff regulators u(k) are considered in our analysis. 
The variation of our final results as u(k) is changed is of order 1% of our determined mass dif- 
ferences, and is included in the quoted error. The deviation as the parameters /, F, C and 6 are 
perturbed is similarly small, and the statistical uncertainty arising from the fit to lattice data is 
smaller still. 

In fact, the dominant uncertainty by an order of magnitude is that that arising from the quoted 
error band on R, the light quark mass ratio. It is clear that better estimates of this value will 
allow our results to be greatly improved in precision, without the need for further lattice data. 
Conversely, a precise determination of the electromagnetic contribution to the n-p mass difference 
could possibly facilitate an improved estimate of R by this method. 

4. Application to QCDSF-UKQCD lattice results 

The method described above can be applied equally to other sets of lattice data for the octet 
baryon masses. In particular, we consider recent 2 + 1 -flavor QCDSF-UKQCD lattice simula- 
tions ffT2l . which follow a significantly different trajectory in the light-strange quark mass plane 
to the PACS-CS simulations. Instead of holding the strange quark mass fixed along the simulation 
trajectory, the QCDSF-UKQCD Collaboration holds the singlet quark mass (m 2 K + ml/2) fixed. 

We use simulation results from this collaboration which lie both along the 'singlet' line and 
along the SU(3) symmetric line [fT2|. Precisely as was done in our analysis of the PACS-CS 
lattice data, we calculate small finite- volume corrections. The lattice spacing a = 0.072(1) fm 
is determined by fixing X N = (1/3)(M N + A/ z + M H ) to the experimental value at physical quark 
masses. This a is somewhat lower than that quoted by the QCDSF-UKQCD Collaboration as 
we account for chiral curvature. The quality of our fit to the isospin-averaged results is clearly 
excellent, with a^ 2 /dof of 0.6 and regulator mass A = 1.0 + 0.1. The fit is shown in Figures [Tj 
[2] and [3} and the values of the octet baryon masses extrapolated to the physical point are given in 
Table|6j These are largely consistent with the physical values. 

Applying the method described in previous sections to the QCDSF-UKQCD lattice data gives 
results for the strong contribution to the octet baryon mass splittings which are consistent with our 
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B 


Mass (GeV) 


Experimental 


N 


0.966(16)(10) 


0.939 


A 


1.112(13)(5) 


1.116 


Z 


1.193(12)(4) 


1.193 




1.307(H)(0) 


1.318 



Table 6: Octet baryon masses based on a chiral extrapolation of the QCDSF-UKQCD data set. The first uncertainty 
quoted is statistical and the second results from the variation of various chiral parameters and the form of the UV 
regulator as described in the text. The experimental masses are shown for comparison. 



fit to the PACS-CS Collaboration simulation results. We find 
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M n - M p = (w/mj )(16.6 + 1.2) MeV (16a) 
M s - - M s + = (co/ml (phy J(4S.9 ± 1.7) MeV (16b) 
M 3 - - M s o = (a/m* )(32.2 ± 1.6) MeV, (16c) 



corresponding to 



M n -M p = 2.4 + 0.3 MeV (17a) 
Afc- - M s+ = 7.0 + 0.9 MeV (17b) 
M s - - M E o = 4.6 + 0.6 MeV, (17c) 

when R is set to the physical value given in Equation [14] 

The QCDSF-UKQCD Collaboration has also recently presented a determination of the strong 
contribution to the baryonic mass splittings, based on these lattice simulations [6]. That calculation 
used a linear and quadratic SU(3) flavor symmetry breaking expansion in the quark masses. As the 
meson and baryon octet expansion coefficients depend only on the average quark mass, provided 
these are kept constant this method allows for an estimation of the baryon mass splittings at the 
physical point using the available 2 + 1 -flavor lattice results. The results of this approach are: 

M n - M p = 3.13 + 0.15 + 0.53 MeV (18a) 
M £ - -M s+ = 8.10 ±0.14 ± 1.35 MeV (18b) 
M 5 - - M s o = 4.98 + 0.10 + 0.84 MeV. (18c) 

The first uncertainty quoted in Equation [T8] is statistical, while the second allows for violations of 
Dashen's theorem. 



5. Discussion 

We have carried out an analysis of i so spin- averaged lattice simulations for octet baryon masses 
using a formal chiral expansion based on broken SU(3) symmetry. Using the resulting expansion 
coefficients one can evaluate the strong contribution to the baryon mass splittings. Our results, 
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based on analyses of PACS-CS and QCDSF-UKQCD lattice data sets, are summarized in Table[7J 
Both calculations yield compatible values, despite significant differences between the two lattice 
simulations, including in particular different lattice sizes, lattice spacings, and different methods of 
determining these spacings. Of course, as emphasized previously, the two simulations also follow 
quite different trajectories in m n - m K space. 

Furthermore, we note that the results of the QCDSF-UKQCD Collaboration analysis of their 
data are are entirely consistent with our values. While the approach taken by the QCDSF-UKQCD 
group makes use of only that lattice data calculated along a trajectory which holds the average 
quark mass constant, we have also included in our fit the data from that collaboration which lies 
away from this line. This contributes to our reduced uncertainties. We also point out that both 
methods require some theoretical input: we input the up-down quark mass ratio R, while the 
the Horsely et al. calculation uses Dashen's theorem (with some uncertainty) to estimate 'pure 
QCD' meson masses at the physical point. The clear consistency between the two independent 
calculations is encouraging. 



A„ u (MeV) 


M„ - M p 


My. - M s + 


M 3 - - M E o 


Chiral (PACS-CS) 


2.9(4) 


7.5(10) 


4.6(6) 


Chiral (QCDSF-UKQCD) 


2.4(3) 


7.0(9) 


4.6(6) 


QCDSF-UKQCD 


3.13(55) 


8.10(136) 


4.98(86) 


Exp. & EM (traditional) 


2.05(30) 


7.91(30) 


5.99(30) 


Exp. & EM (WLCM) 


2.6(6) 







Table 7: Up-down quark mass contribution to octet baryon mass splittings. Lines 1 and 2 show the results of our 
chiral extrapolation of PACS-CS and QCDSF-UKQCD lattice data respectively. Line 3 shows the QCDSF-UKQCD 
Collaboration analysis of their data as described in the previous section, while lines 4 and 5 give estimates deduced 
from the total mass splittings and electromagnetic contributions, as determined by Gasser and Leutwyler (traditional) 
or Walker-Loud et al. (WLCM). 

While our results and those of the QCDSF-UKQCD Collaboration are consistent with both the 
traditional and Walker-Loud et al. (WLCM) determinations of the strong contribution from the 
electromagnetic component (see Table[7]), we note that there is a small discrepancy between our re- 
sults and the traditional value for the cascade baryon mass splitting. Clearly, determinations of all 
octet baryon electromagnetic mass splittings using the WLCM analysis would be of considerable 
interest. 

While more lattice data for isospin-averaged octet baryon masses, on larger lattice volumes, 
would allow the uncertainties of our calculation to be somewhat reduced, we emphasize that the 
dominant uncertainty in our calculation arises from the up-down quark mass ratio. A more precise 
value of m u lmd could reduce the uncertainty of our determination by as much as an order of 
magnitude. Conversely, direct lattice determinations of the electromagnetic contributions to the 
mass splittings, with the analysis presented here, may act to significantly improve the value of 
m u /m d . 
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Figure 1: Fit to the QCDSF-UKQCD baryon octet data. Error bands shown are purely statistical, and incorporate 
correlated uncertainties between all fit parameters. Note that the data shown has been corrected for finite volume 
and the simulation strange quark mass, which was somewhat different from the physical value at each point. The red 
circles and green squares lie on the singlet trajectory and the SU(3) symmetric line respectively, and the blue stars 
denote the physical points. 
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Figure 2: Fit to the QCDSF-UKQCD baryon octet 
data, shown along the SU(3) symmetric line. Error 
bands are as in Figure [T] 



Figure 3: Fit to the QCDSF-UKQCD baryon octet 
data, shown along the singlet trajectory. Error bands 
are as in Figure [T| 
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